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PREFACE 



This Technical Report is based on the master's thesis of Tom R. Houston. 
Members of the thesis committee were Julian C. Stanley, Chairman; Frank 
B. Baker; Henry B. Mann; and J. Marshall Osborn. 

Extending knowledge about, and improving educational practices related to, 
cognitive learning in children and youth is the primary goal of the Wisconsin 
R & D Center for Cognitive Learning. The Laboratory of Experimental Design, 
part of the technical section of theR & D Center, provides valuable assistance 
to project directors in the design of experiments and also in the analysis of 
data. Further, the staff of the LED are charged with extending knowledge about 
experimental designs, scaling procedures, data analysis and the like. 

This Technical Report is the fifth in a series describing new developments 
in the methodological area. It deals with Latin squares as a control for pro- 
gressive and adjancy effects in experimental designs. The history of Latin 
squares is reviewed, and several algorithms for the construction of I^atin and 
Greco-Latin squares are proposed. Results are of particular application to 
rating studies, and to designs requiring the "rotation" of teachers in class- 
room research. 
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INTRODUCTION 



Latin. sQuareSi variously called, "magic squares" and "Eulerian squares »" are of considerable 
antiquity. The epigraph of this paper gives an example from astrology that may go back to Hellen- 
istic times, and the extension of this use into ritual magic is well knovn; e.g., in Marlowe's (10) 
play The Tragedy of Doctor Faustus (1604), reference is made to 

. . . J ehovah' s name 

Forward and backward anagrammatized (3. 9, 10) 

probably alluding to a Latin square, k = 4, constructed by (19. 2) below, whose types were the 
Tetragrammaton. It is reported (10) that this square summons the Devil, for which reason it will 
not be exhibited here. 

In recent times Latin squares have been found to have further practical application in experi- 
mental design. Fisher ( 6 ) places the use of the Knut Vik square (a permutation by rows of the 
square described in (6. 2) below) in agricultural experimentation as early as 1872. It is con- 
venient to categorize applications of Latin squares in experimental design under three headings: 

I* To increase the number of experimental units . A field subjected to k treatments can be 
divided by restricted randc j lization into k^ plots by a Latin square design, giving a k-fold increase 
in experimental units (and distributes soil variability) for the same acreage. 

II. For fractional factorial designs . Latin and hyper-Latin squares are sometimes used to 
sample possible treatment combinations within k”*' fractional designs. Since such designs involve 
k*' aliases for every effect by the confounding of interactions, their use is questionable unless 

there is assurance, rarely available in the behavioral sciences, that the interactions are close to 
zero. 

Ill* For counterbalancing repeated measures designs . When each experimental subject re- 
ceives k successive treatments, it is often feasible and desirable to equalize progressive effects 

(warm-up, fatigue, practice, etc.) by Latin squares. These insure that every treatment occurs, 

1 

on the average, in the 2 (^+l)th ordinal position, and in every position with equal frequency. 

Note that the last application is a case of the second, where ordinal position is considered as 
a non-interactive factor with k levels. A difficulty can arise in such cases, as the following 
hypothetical example illustrates: 




IX 



Suppose that a study of diabolism were to be conducted, to determine the relative effective- 
ness of four methods of exorcism; Bell (T^), Book (T2)» Candle (T3), and Team Exorcism (T4). 
The dependent variable could be the response to the question "Do you renounce the Devil?", di- 
chotomously scored (affirmative = 1 , no response or negative = 0 ). Suppose further that the in- 
vestigator wished to apply all treatments to each S (Ss being 4 N demoniacs), employing a Latin 

square design of the form of ( 9 . 1 ) below, so that non-ultimate T^ was always followed by T2» T2 
by T3, T3 by T4, and T4 by T^. Then if T2 elicited a significantly higher mean response than 
Ti, there would still remain several plausible alternative hypotheses to the interpretation that 
T2 was more potent than Tj^. Perhaps T4, which preceded Tj^ in 3 N cases, left Ss unable to re- 
spond until after T2; or perhaps T2 is comparatively ineffective unless an immediately prior Tj^ 
has rendered Ss susceptible to its effects. 

Such "residual effects" are a general problem in repeated measures designs which an injudi- 
cious choice of Latin squares can aggravate, as in the above example. Randomization (i. e. , ran- 
dom linear permutations by row and column of a randomly selected square) is often suggested as 
a remedy, but in dealing with small numbers there is considerable risk to the investigator of sub- 
stituting bad luck for bad judgment and obtaining a design with unfortunately biased residuals con- 
taminating main effects. Complete sets of permutations avoid this, but there may not be enough 
Ss for this to be possible. 

As Grant (8) and others have observed, such effect cannot be entirely controlled by Latin 
squares. It is often worthwhile, however, to control for the most potent of these, the immedi- 
ate sequential effects, by constructing Latin squares in such a manner that every treatment is im- 
mediately followed by every other an equal number of times. Existing procedures for accom- 
plishing this will be reviewed in this paper, and a general solution offe^’ed for the previously un- 
solved case of k odd, where k is the order of the Latin square. This solution also provides for 
the counterbalancing within (but not across) alphabets of Greco-Latin squares for k odd. It 
should be noted that the given procedure ( 23 . 1 ) counterbalances only within rows; the extension to 
within-column counterbalancing is trivial, but seems to be without application in the behavioral 
sciences. In the diabolism example given above, incidentally, Faustus' square ( 19 . 1 ) will coun- 
terbalance immediate sequential effects. 
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PART I 

ORTHOGONAL SQUARES 



Throughout thij discussion, k will be the symbol for the order of a Latin square (defined 
below). For some natural number c, k' = 2c + 1, and k" = 2c. The notation e-. (sometimes 
with a comma in the subscript) will represent a cell in rectangular array E, at the intersec- 
tion of the ith row and the jth column, i and j ranging from 0 to (k- 1 ). E- represents the 
type occupying cell e^j, and E^^ the type in q* Thus ©2,5 “ ^7 signifies that the same type 
is in cell 5 0* 

Definition 1 ; Rectangular array E is a Latin square if: 

( 1 . 1 ) E.j Epj for all i / p, j / q. 

(1.2) There are k = i(max)+l = j(max)+l t%pes. 

(1.3) Exactly one type occupies each cell. 



The number k is called the order of square E. An example of a Latin square of order k = 
5 is: 

(1.4) 1 3 4 0 2 

4 2 0 3 1 

2 1 3 4 0 = (1.4)E. 

0 4 2 1 3 

3 0 12 4 



A Latin square will be called a standard form (SFL) if the types in its first row and col- 
umn are in natural order. Any Latin square can be put into the form of a SFL by pre- and 
post-multiplication with appropriate permutation matrices. For example, the square given 
above may be standardized by the operation: 



( 2 . 1 ) 



0 0 0 1 
10 0 0 
0 0 10 
0 0 0 0 
0 10 0 



0 

0 

0* (1.4)E 
1 
0 



10 0 0 

0 0 0 0 

• 0010 

0 10 0 

0 0 0 1 



0 0 1 2 3 4 

1 1 0 4 2 3 
0=24301 
0 3 2 1 4 0 
0 4 3 0 1 2 . 



Transformations such as the above, changing the order of the rows and/or columns of a 
Latin square, will be called linear permutations. Call all the SFLs obtainable from a given 
Latin square a linear set (a term not in general use). A number of higher-order groupings 
are available for the classification of Latin squares [see Norton (11)* Fisher and Yates (tD. 
including transformation sets, species, families, domains, only the first of Which will be 



discussed here. An interesting feature of all classificatory systems, however, is that as of 
the present, the cardinality of the set of all SFL.s, linear sets, species, etc. is unknown for 
any k > 7. 

Two Latin squares of the same order, R and G, are said to be orthogonal to each other if , 
when juxtaposed by corresponding elements to form a k order array of double-entry cells 
(R,G), every paired combination of types exactly once. (Latin squares are 

so named, incidentally, because Roman letters are often used to represent the several types 
of such squares. Similarly, squares of the form of (R,G)are called "Greco-Latin squares" 
because Greek letters are often used to represent the types of one of the orthogonal squares. ) 
An example of a Greco-Latin square for k = 3 is; 



(3.1) 



00 


12 


21 


11 


20 


02 


22 


01 


10. 



Orthogonal squares cannot be constructed for k = 2 or 6. Gunther (14) cites an early proof 
of the latter (1842) by Clausen by exhaustive examination of 6x6 squares. Akar (12), also 
cited in Norton's extensive review of the literature before 1939. gave aproof in 1895 thator- 
thogonal squares existed for any k' odd, doubtless along lines similar to those followed in 
Algorithms I or II below. The even numbers proved more refractor , Bose (3) and Stevens 
(17) independently and concurrently proved the existence of orthogonal squares fork = 2*' and 
k = p'*' in general, for p' a prime and n a natural number greater than one. Mann (9) pro- 
vided for the construction of non-cyclic (see section 6) orthogonal squares for a lax'ge class 
of values of k; but not until 1959 did Bose and Shrikhande (4) disprove Suler's (13) conjecture 
that no orthogonal squares exist for k = 4t+2, for t an integer. 

For k = 5 there exist 56 SFLs. These fall into sixteen linear sets, ten of five SFLs and 
six of one. Appendix A lists these SFLs in lexicographic order, indicating the transforma- 
tion set to which they belong. To illustrate the concept of a transformation set, define Q(E), 
called a recoding of a Latin square E, as the square constructed by substituting throughout 
E Qq for Eq, Qi for Ep. . . , for Ej^^p where Qq.Qi,. . . any permutation of 

the first k non-negative integers. Call $(w) any linear permutation of a SFL (for k = 5) des- 
ignated in Appendix A as a member of linear set (w). Then the notation 

(4.1) SFL(l) 02413 = $(6) 

indicates that any SFL from set (1), when recoded by the permutation 02413, gives a square 
whose SFLs are in set (6). The reader can easily verify that: 



SFL(l) 


13024 = $(7) 


SFL(l) 


01234 = $(1) 


SFL(l) 


24130 = $(8) 


SFL(l) 


12340 = $(2) 


SFL(l) 


30:*41 = $(9) 


SFL(l) 


23401 = $(3) 


SFL(i) 


41.-i02 = $(10) 


SFL(l) 


34012 = $(4) 






SFL(l) 


40123 = $(5). 



For any integer m 0, if x = am + r, for integer a and 0 - r < m, x is said to be congru- 
ent to r modulo m, written x^ r(mod. m), and x = r (reduced modulo m). Define V as the 
union of linear sets (1) through (5), and V' as the union of linear sets (6) through (10), v and 
v' being members of these sets. The reader can verify that 

(4.3) SFL(v) (0+y, 1+y, 2+y, 3+y, 4+y) (mod. 5) = $(v) 

SFL(v') (0+y, 1+y, 2+y, 3+y, 4+y) (mod. 5) = $(v') 

SFL(v) (0+y, 2+y, 4+y, 1+y, 3+y) (mod. 5) = $(v') 

SFL(v') (0+y, 2+y, 4+y, 1+y, 3+y) (mod. 5) = $(v), 

and analogous statements regarding V and V' can be made for similar sets of permutations, 
within each of which relationships such as occur in (4. 1,2) exist. Since by definition, Latin 
squares within the same linear set may be permuted into each other, by appropriate recod- 
ings any square in the union of V and V' can be transformed into any other. A transforma- 
tion set is defined as a set of all Latin squares that by recodings and linear permutations can 
be changed into each other. 



It is not difficult to show that the six 5x5 SFLs in Appendix A not designated as members 
of linear sets (I)-(IO) are in fact not members of the transformation set V union V'. Con- 
sider the squares 



(5.1) 



0 12 3 4 

1 0 3 4 2 

2 3 4 1 0 = T, 

3 4 0 2 1 

4 2 10 3 



0 12 3 4 

1 4 0 2 3 

2 0 3 4 1 = C, 

3 2 4 1 0 

4 3 10 2 



where T is a member of the union of linear sets (I)-(IO), while G is not. It is claimed that 
no recoding or linear permutation of T will give C. 



Note that in T there exist "rectangles" of the form; 



(5.2) 




^i,j+n 



^i+m, j 



t 



i+m, j+n 



such that and for n,m 0. An example of such a rectangle 

in T is where ^0,3’ ^0 1 ~ ^2 3’ i~0, j = l, m = n = 2. Clearly in any recoding 

Q(T) of T, q(t)Q j = q(t)^ 2 ^(t)^ ^ ~ 3* Similarly the rectangles persist over lin- 

ear permutations, which merely add constants (modulo k) to i,j,m, and n. The square C has 
no such rectangles, so C is outside the transformation set of V union V'. By Algorithm I 
below, it can be demonstrated that the six 5x5 Latin squares outside of V union V are all 
members of a second transformation set. 



For k - 5, the SFLs having no such rectangles (intercalates is another name forsuchcon- 



figurations) can be recoded and linearly permuted into the SFL; 
(6. 1) e^j = (i + j) (mod. k). 



Such squares are called cyclic Latin squares (CLS), and can be constructed for any natural 
number k. For any columnwise linear permutation of a CLS, it immediately follows from 
this definition that 

(6. 2) if eij = and e^ j+p = E^+q. 

then 

where all subscripts (as throughout this paper) are reduced modulo k. 

Note that the orthogonal squares in (3. 1) are both columnwise permutations of theCLSfor 
k = 3. Greco-Latin squares of any order k' can always be obtained fromCLSsby the follow- 
ing construction: 

Algorithm I; Prepare an empty k'xk' matrix, k' = 2c+l. Index the first c columns aj,a 2 ».. 
. , a^ from left to right and the last c columns * '^1 right. Index the 

central column Assign k values Vj to these column indices so that these values form a 

permutation of the first k non-negative integers with the restriction 

(7. 1) (aj - bj) ^ (aj^ - b^^) ^ (bj - aj) for all j ^ n, 

all values reduced modulo k. In every column j let 

(7.2) e-. = (i - V.) (mod. k). 

J 

Call the square thus constructed A. Construct square M: 

(7.3) 

where a^^^ j = e^^j in (7. 2). It is claimed that A and M are orthogonal Latin squares, and that 
any M constructed from Q(A) is orthogonal to Q(A). 

Example of Algorithm I for k' = 7: Let Vj|^ be the permutation 1564302, which satisfies (7. 1). 
Constructing A by (7. 2), and recoding by the permutation (randomly chosen) 3401265 gives 
Q(A) = 

(8. 1) 3 1 0 2 6 4 5 

4 2 1 6 5 0 3 

0 6 2 5 3 1 4 

1 5 6 3 4 2 0- 

2 3 5 4 0 6 1 

6 4 3 0 1 5 2 

5 0 4 1 2 3 6 

The reader can verify that in (8. 1) Q(A) Q(M) 

That A will always be a Latin square (from which it follows that M, Q(A), and Q(M) are 
also) follows immediately from (7. 2) and from the specification that the values ofvj^be a per- 
mutation of the first k non-negative numbers; A will always be a columnwise linear permu- 
tation of the CLS, as for example; 



0000001 

0000010 

0000100 

0001000 = Q(M). 
0010000 
0100000 
1000000 



(9.1) 



0 k-1 k-2 ... 1 

1 0 k-1 ... 2 

2 1 0 ... 3 



k-1 



k-2 k-3 



0 . 



The same construction (9. 1) serves as a proof that (7. 1) always has at least one solution for 
any k'. In (9. 1), forms the permutation 0, 1 , 2, ..., k-1 , so that 

(9. 2) (a. - b.) = (2j - 1) (reduced modulo k) 

J J 

(9. 3) (bj - aj) = (1 - 2j) (reduced modulo k). 

Since k' = 2c + 1, and j - c, the right-hand term in (9. 2) is an increasing function of j, and 
in (9. 3) a decreasing function of j as j varies from 1 to c. And since these are of opposite 
parity (the first odd, the second even) it follows that (7. 1) is satisfied by (9. 1). 



To prove that A J_ M, for all A constructed by Algorithm I, it shall be shown that (Q(A), 
Q(M) is always a Greco-Latin square. In the first c columns of (Q(A),Q(M)) entries are of 
the form: 



(10.1) (A^.An + aj - bj), 

the subscripts referring to the first column of Q(A). In the last c columns the corresponding 
expression is: 

(10.2) (Aj^,A^ + bj - aj). 

In the median column the expression is: 

(10.3) (An, An). 

(7. 1), for a given n, ( 10. 1 , 2, 3) assume k different values in the k columns of (Q( A), Q(M) ). 

But n has k different values; therefore (Q(A), Q(M)) has k^ different entries, which proves that 

it is a Greco-Latin square. 

Note that for k = k", m- a^^j by (7. 3), so in (Q(A), Q(M)) entries of the form of (10. 3) can- 
not occur. This proves that the algorithm fails for k even. Euler (13) has in fact shown that 
for k” no Latin square orthogonal to a CLS (or any permutation or recoding thereof) exists. 

Any assignment of Vj^ values satisfying (7.2,3) can be represented as a circle having k' 
equally spaced points on its perimeter, 2c of these being connected by c chords of c different 
lengths. These points are identified with the first k non-negative integers, numbered con- 
secutively clockwise or counterclockwise from any point. For k' = 9, 



is such a figure. Identify each chord with a different number j, so that the values of j form 
a permutation of the first c natural numbers, and assign to one endpoint of every chord the 
symbol a, and b to the other endpoint. Thus each chordal endpoint has a unique name aj or 
bj. Assign to the point lying on no chord the name By construction (7. l)is satisfied by 

a permutation of the first k non-negative integers. For k' < 50, this is in fact a convenient 
method of solution, alternative to (9* 1). In experimental designs, (A,M) can be randomly 
permuted and lineraly recoded, without affecting orthogonality. 

12. For k' = ab. Algorithm I can be iteratively used to construct squares like: 



A 


B 


C 


D 


E 


F 


G 


H 


I 








B 


C 


A 


E 


F 


D 


H 


I 


G 


(R) 


(S) 


(T) 


C 


A 


B 


F 


D 


E 


I 


G 


H 








D 


E 


F 


G 


H 


I 


A 


B 


C 








E 


F 


D 


H 


I 


G 


B 


C 


A : 


= (S) 


(T) 


(R) 


F 


D 


E 


I 


G 


H 


C 


A 


B 








G 


H 


I 


A 


B 


C 


D 


E 


F 








H 


I 


G 


B 


C 


A 


E 


F 


D 


(T) 


(R) 


(S). 


I 


G 


H 


C 


A 


B 


F 


D 


E 









Call M constructed by (7. 3) the mirror image of A, for any Latin square A. The reader can 
verify that the square given in (12. 1) is orthogonal to its mirror image. For the case of dou- 
ble use of Algorithm I, if k' = gh, it will be possible to partition k types into g sets of h; call 
these sets A, B, . . . , G. Using the types in each set, construct by Algorithm I g Latin 
squares A, B, . . . ,G of order h, using the same values for Vj^ in each square. 

Now treating A,B,...,G as elements in a super-column construct by Algorithm I a Latin 
super-square letting k' = g and c = ■j(g-l). Note that these constructions are always possi- 
ble, since k' = gh = 2c + l implies g,h ^ l(mod. 2). The resulting super-square is not in gen- 
eral a CLS permutation, and can easily be changed into k' order squares of different trans- 
formation sets having orthogonal mirror images by the judicious rearrangement of types with- 
in one or more of the h-order squares. 

13. To digress at some length, let it be remarked that by a construction quite similar to Al- 
gorithm I orthogonal Latin squares can be obtained for any k' odd, one of which is not nec- 
essarily in the same transformation set as the CLS. Definition : A diagonal of square array 
E is any set of cells in E of the form: 

(13.1) e.^Q, e.^2,2> •••> ®i+k-l,k-l* 

In this section, unless otherwise indicated, all values are reduced modulo k. From section 
1 it follows that a Latin square has k diagonals of k cells each. Note that any CLS put into 
the form of (9. 1) has only one type on each diagonal, while its mirror image has k types on 
each diagonal. It follows at once that (9* 1) is orthogonal to its mirror image. 

Algorithm II ; Prepare an empty k'xk' matrix, and obtain some permutationDQ, , . . . 
of the first k non-negative integers, such that; 
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(13.2) (D^ - Dq) (i - q) for all i q. 



Write this permutation into the first column of matrix D. In every column j, let: 
(13.3) dj. = (D._.) + (j). 



This implies that: 




If C is a CLiS in the form of (9. 1), it is claimed D J_ C. 

Since orthogonality follows from the construction, it shall only be shown that D is a Latin 
square, or specifically, that (1.1) holds, or that 

(13. 5) Djq i Djj + Dpj for all i p, j q. 

Suppose Dj^q = D- = Dg, for some z,i, with j q. Then 

(13.6) °i-j,j-j = (13.4), and 

l^i-q,q-q " (^^^ therefore 

i^i-j - i^i-q = (q - j) 

which contradicts (13.2). Therefore ^ D-. Suppose D- = Dpj = for some z,j with 
i p. Then by (13. 4) 

(13.7) Dp-j.j-j = (°z) - (j)j by (13.6) 

which contradicts the definition of Dq,Dj,... Then (13.5) is proved, andD is a Latin 

square. 

For k' = 7, Algorithm II gives the square: 

(13.8) 2150364 0654321 



It can be proven that D in (13. 8) is not in the CLS transformation set. Note that i = 
^4,4' ^2,4 ~ ^4, 1' ^ rectangle of the form of (5. 2). If such a rectangle could exist in a 

a CLS, by (6. 1, 2) E^^+q = and Ej^+^^+q = Ej^ for some q 0, E a CLS. This implies 

thatm = q,and that m+q = k, or that 2q = k. Since q is an integer and k odd, this is impossible. 

It is stated without proof that (13. 2) is satisfied for any k' by: 

(13.9) l^i “ (yl “ 2 ) (mod. k) 

for any z and any y > 1 where (y,k) = 1. This rule always gives a D in the CLS transforma- 
tion set. In (13. 8) the permutation 2516430 gave a non-CLS D, but no general method of ob- 



5 3 2 6 1 4 0 

1 6 4 3 0 2 5 



1 0 6 5 4 3 2 

2 1 0 6 5 4 3 

3 2 1 0 6 5 4 

4 3 2 1 0 6 5 

5 4 3 2 1 0 6 

6543210 . 



6 2 0 5 4 1 3 = DJ_C 



4 0 3 1 6 5 2 
3 5 1 4 2 0 6 
0 4 6 2 5 3 1 
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taining such permutations is known to this writer. Appendix B lists several non-CLS 
permutations. 

It ia natural to ask whether Algorithm II gives orthogonal squares for k even; herein might 
lie a simple counterexample to Euler's conjecture (see section 3). But like Algorithm I, Al- 
gorithm II can be shown to fail for all k = h". Suppose Dq,Dj^, . , . were a permutation 

of the first k non-negative integers satisfying (13. 2). Define = i - D- for i = 0, 1, 2, . . . , 

k"-l. It follows from (13. 3) that bj^ = bq for i q. Since bj^ is reduced modulo k", this im- 
plies that: 

(13.10) i=k"-l i=k"-l 

“ Y/ ^ ~ ^c(mod. k"). 

i = 0 i = 0 

But by the definition of bj and by the fact that the Dj^'s are a permutation of the first k non- 
negative integers, 

(13.11) k"-l k"-l k"-l k"-l 

Z) ^i “ Z) ^ ” Z; ^i “ 

0 0 0 0 

which contradicts (13. 10), so Algorithm II fails for k". 

Despite this result, Euler's conjecture can be disproved by a simple combination of tech- 
niques used in Algorithms I and II. For example, Parker's (15) 10x10 Greco-Latin square, 
constructed by an extension of the methods of Bose and Shrikhande (4), is composed of a 7x7 
square with ten types from each alphabet orthogonal by diagonals; a 7x3 and a 3x7 rectangle 
with cyclically orthogonal rows and columns respectively, each with seven types from each 
alphabet; and a 3x3 square of three types from each alphabet, orthogonal by any method. Ap- 
pendix C gives an 18x18 square analogously constructed, and the outline of a general method 
for obtaining anti-Eulerian squares. See Barra and Guerin (2) for k" ^ 6(mod. 12). 



PART II 

SEQUENTIAL COUNTERBALANCING 



14. Every row i of a Latin square has k-1 pairs of consecutive elements of the form Eij 
Ei,j+1- 

(14.1) 4 10 2 3 

from a square where k = 5, has four pairs: 4 1; 1 0; 0 2; and 2 3. Since there are k rows, 
there are k(k-l) pairs in E of order k. Defining E^ E^ ^ it can be seen that k(k-l) is 

the number of possible pairs of types in E. Definition ; A Latin square containing all possi- 
ble consecutive pairs in its rows is said to be counterbalanced for immediate sequential ef- 
fects (abbreviated SQCB). 

15. Theorem 1 ; A square in the CLS transformation set is SQCB if and only if the set of all val- 

ues d(j) is a permutation of the first k-1 natural numbers, where d(j) = (s - r)(mod.k), E^^j 

= Ej., E^^ - Eg. 

From (6. 1) it follows that d(j) is constant for any i, despite recodings or linear permutations 
by rows or columns (although these introduce new constant values). Suppose E were a CLS 
such that d(j) = d(j>'-‘) for j j*. Then E^^^ (= Ej^) is followed by But in 

some row i i, E^,.^ = E^^^ is followed by E^,;, + i = E^+jj(j). So E^^^ 

occurs twice in E. Then if E were SQCB the remaining k(k-l)-l pairs of types would occur 
in only k(k-l)-2 pairs of cells, which is impossible. Thus the theorem is necessary. 

Proof of sufficiency follows from the theorem's requirement that type Ej^ be followed by 
types Ejyj^p Ej^^ 2* • • • * ^m+k-1 every E^^^ occurring in the firstk-1 columns of E. Re- 
duced modulo k, this gives k-1 different types following E^« Therefore all k types in E are 
followed by k-1 types exactly once. 

Corrolary ; No CLS-transformation set square is SQCB for k = k'. 

In such a square, if it existed, by Theorem 1: 

(15.1) j=k-2 i=k-l 

^ d(j) = Z) ^ ” ’I’Mk-l) = kc ^ 0(mod. k). 
j = 0 i = i 

But then by (6. 1,2), in such a square where E^^ q ” ^m* then by the definition of d(j), 

= = Ejjj, which contradicts (1, 1). 
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Bradley (5) gives a method of constructing SQCB permutations of CLSs for k” that satis- 
fies the restriction of Theorem 1. He specifies that d(j) = (-l)J(j), by letting the first row be 
of the form: 



(16«1) ^ 0 * ^k“ 1 *^ 1 * ^k- 2 * ^ 2 * * * * * ^c “1 > ^c +1 * ^c * 

and setting the types equal to the subscripts in (16. 1), 

(16.2) =(£(,_.)+ ( 1 ) 

describes the rest of the cells. Since for j' ^ l(mod. 2 ), d(j') is a monotonically decreasing 
function of j, and for j" ^ 0(mod. 2 ), d(j”) is monotonically increasing, then Bradley's d(j) 
forms a permutation of the first k -1 natural numbers, d(j") and d(j') being of opposite parity. 

For k = p'-l, another SQCB CLS is available: 

(17. 1) 123456 123456 

246135 362514 

362514 246135 

4 1 5 263 = S;S»!‘ = 65 432 1 

531642 415263 

654321 53164 2; 

S is an example for p' = 7 (p' must be a prime). S’!* is a linear permutation of S by rows, by 
which it can be seen that S is a recoded and linearly permuted CLS. In S, the SFL, such 
squares are symmetrical around both diagonals andhave the peculiarity that if ©ij ” ej^ = 
and E^+rj+q = Ej^, then ^^ 4 .^ where all subscripts are reduced modulo p'. 

S is constructed by the rule 



(17.2) 






IJ “"(j*i)’ 

again reducing by modulo p', where i and j range from 1 to k, and are in fact simply modu- 
lar multiplication tables for p'. Alimena ( 1 ) was apparently unaware of this when he described 
his procedure for constructing them (which might also serve as a masterplan for a spider- 
web). Alimena asserts that squares of the formof (17. 2 ) are superior for experimental work, 
in terms of sequential properties, to those described by Bradley (I 6 . 1, 2 ). 



Define the separation between two types in a row of a Latin square as the absolute differ- 
ence between the column indices of the cells in which they appear.. The mean separation is 
the summed separations for two letters divided by k. (Since this is to bean ad hoc means of 
comparing (17. 2) and (16. 1,2), the directional insensitivity of thiese definitions is unimpor- 
tant; by both constructions every row is the mirror image of some other row. ) For k = 10, 
mean separations by types are tabled in (18. 1). Entries above the main diagonal have ref- 
erence to ( 16 . 1 , 2 ); those below, to (17.2), or rather to Alimena's construction, which has 
identical separation values. Note that for an ideally counterbalanced sq\xare, mean separa- 
tions would be constant across all types. 



(18.1) 





Eq 


El 


Ez 


E 3 


E 4 


E 5 


E 6 


E 7 


Eg 


^9 


Eo 


- 


1.8 


3. 2 


4. 2 


4.8 


5. 0 


4.8 


4. 2 


3. 2 


1.8 


El 


3. 0 


- 


1.8 


3. 2 


4. 2 


4.8 


5. 0 


4.8 


4. 2 


3. 2 


Ez 


3. 2 


3.8 


- 


1.8 


3. 2 


4. 2 


4.8 


5. 0 


4.8 


4.2 


E 3 


3. 2 


3. 0 


4. 0 


- 


1.8 


3. 2 


4. 2 


4.8 


5. 0 


4. 8 


E 4 


4. 0 


3.8 


4. 0 


3. 2 


- 


1.8 


3.2 


4. 2 


4. 8 


5. 0 


E 5 


3. 0 


3. 2 


3. 0 


3.8 


5. 0 


- 


1.8 


3. 2 


4.2 


4.8 


E6 


3. 8 


4. 0 


3. 0 


5.0 


3. 8 


3. 2 


- 


1.8 


3.2 


4. 2 


E 7 


3. 8 


3.2 


5. 0 


3. 0 


3. 0 


4. 0 


4. 0 


- 


1.8 


3.2 


Eg 


4. 0 


5. 0 


3. 2 


4.0 


3. 2 


3.8 


3. 0 


3.8 


- 


1.8 


E 9 


5. 0 


4. 0 


3. 8 


3. 8 


3. 0 


4. 0 


3.2 


3. 2 


3. 0 


- 



The average absolute deviation from the grand mean is 0.5 for (17.2), 1.0 for (16.1,2). 
For a 10x10 square, the greatest possible mean separation is 5.0, the least 1.8; (16.1,2) 
achieves this range, and (17.2) the upper limit. Since it is generally less detrimental that 
successive treatments tend to be remote than that they tend to follow one another closely, 
(18. 1) supports the view that modular tables give "better” squares than (16. 1,2); but the dif- 
ference may not be of practical importance, and it must ' e remembered that (17. 2) is only 
applicable when k = p'=l. 



19. It can be shown that (17. 2) will not generalize to 



(19.1) ®ij = E(j-i)» 

with the subscript reduced modulo k+1 p'. If E were a square so constructed, then k+1 = ab 

for some a, 1 < a - k^, and in some cell where i, j range from 1 to k, ej^^j = 

= Eq. But i* 1 0(mod. (k+1)) in the range of i, so Eq E^, l so E is not a Latin square. 



20. By the corrolary to Theorem 1, no SQCB Latin square exists in a CLS transformation set 
for k'. No SQCB Latin square in any transformation set is known fork', and it is conjectured 
that none exists. 



It will be shown, however, that for k' it is always possible to construct two Latin squares 
E and F such that each pair of k(k-l) possible pairs of types occurs exactly twice within the 
rows of E and F. Such squares will be said to be complementary. 

Algorithm III; For k' = 2c + 1, squares E and F where 

(20.1) ey = (-l)j Ei + 1)/|]+ i 

can be constructed, where Qc] represents the integral part of x, such that E is complemen- 
tary and orthogonal to F. 
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It is evident that each column of E and F conforms to (6.2). ThenEandFare cyclic Latin 
squares if for constant i, ej^j and fj[^j form permutations of the first knon-negative integers 
as j ranges from 0 to k-1. Note that for j” even, in (20. l)is a monotonically increas- 

ing function of j"; and for j' odd, j' is a monotonically decreasirg function of j'. But 

(20.2) c - ei,j"(niax) “ ®i,j'(min) = c + i» 

an analagous statement holds for j* Then by monotonicity, E and F are Latin squares. 

To show E J_ F, set 

jj + c + 1 for j even; 

(20.3) b(j) = < 

Jc - j for j odd. 

Then Ej^j - F j ^b(j) (mod. k). Now if b(2t) ^b(2t'+l) (mod. k), then 2t + c + 1 ^ c + 2t' - 1 
(mod. k), so 2t + 2t' + 2 0(mod,k). But since 2t + 2t' + 2 < 2k, v/e must have 2t + 2t' + 2 = 

k, which is impossible, since k is odd. Since (20. 3) gives two monotonic functions, this 
proves that E J_ F. 

Furthermore, for d(j) in E, as defined in section 15 above, d(j) assumes the values 2,4, 
. . . ,k-l twice. Then every pair of types Ej^^ j is of the form 

(20.4) Ej^ Ej^^2y some natural number y. 

It follows that kc ordered pairs of types occur in E. Similarly, F contains kc pairs of the 
form 



(20.5) Ejj^ Ej^4.2y_i 



since d(j) assumes c different odd values twice in any row. Sincr ;20. 4) f (20.5), it follows 
that kc + kc = k(k-i) different pairs occur twice in E and F, so E and F are complementary. 



for k = 7; 



( 21 . 1 ) 0 6 
1 0 
2 1 

3 2 

4 3 

5 4 

6 5 



3 

4 

5 

6 

0 

2 1 
0 4 13 2 



3 4 

4 5 

5 6 

E; 6 0 

0 1 
1 2 
2 3 



2 5 16 0 

3 6 2 0 1 

4 0 3 1 2 

5 14 2 3 

6 2 5 3 4 

0 3 6 4 5 

1 4 0 5 6 



Example of Algorithm III 

15 2 4 

2 6 3 5 

3 0 4 6 

4 15 0 

5 2 6 1 

6 3 0 



F. 



The extension of this construction to within-column counterbalancing is accomplished by 
permuting the last k-1 rows until the first column is the transpose of the first row in both 
squares. The construction will also give pairs of Greco-Latin squares counterbalanced for 
sequential effects within (but not between) alphabets. 



APPENDIX A 
THE 5x5 SQUARES 



The k = 5 SFLs in Lexicographic Order ; 

First rows and columns are deleted; linear sets indicated. 



(1) 






(10) 






(5) 




(9) 






(9) 






(5) 






(1) 






0 3 


4 


2 


0 3 


4 


2 


0 3 4 


2 


0 3 


4 


2 


0 4 


2 


3 


0 4 


2 


3 


0 4 


2 


3 


3 4 


0 


1 


3 4 


0 


1 


4 0 1 


3 


4 1 


0 


3 


3 0 


4 


1 


3 1 


4 


0 


4 3 


0 


1 


4 1 


2 


0 


4 0 


2 


1 


2 4 0 


1 


2 4 


1 


0 


4 1 


0 


2 


4 0 


1 


2 


2 1 


4 


0 


2 0 


1 


3 


2 1 


0 


3 


3 1 2 


0 


3 0 


2 


1 


2 3 


1 


0 


2 3 


0 


J 


3 0 


1 


2 


(10) 






(4) 






(3) 




(8) 






(7) 






(8) 






(7) 






0 4 


2 


3 


2 0 


4 


3 


2 0 4 


3 


2 0 


4 


3 


2 0 


4 


3 


2 3 


4 


0 


2 3 


4 


0 


4 3 


1 


0 


3 4 


0 


1 


3 4 0 


1 


4 3 


0 


1 


4 3 


1 


0 


0 4 


1 


3 


0 4 


1 


3 


2 0 


4 


1 


4 1 


2 


0 


4 1 0 


2 


0 4 


1 


2 


0 4 


2 


1 


4 0 


2 


1 


4 1 


0 


2 


3 1 


0 


2 


0 3 


1 


2 


0 3 2 


1 


3 1 


2 


0 


3 1 


0 


2 


3 1 


0 


2 


3 0 


2 


1 


CLSa 




(9) 






(6) 




(10) 






(4) 






(3) 






(5) 






2 3 


4 


0 


2 3 


4 


0 


2 3 4 


0 


2 3 


4 


0 


2 4 


0 


3 


2 4 


0 


3 


2 4 


0 


3 


3 4 


0 


1 


4 0 


1 


3 


4 1 0 


3 


4 1 


0 


3 


0 3 


4 


1 


0 3 


4 


1 


3 0 


4 


1 


4 0 


1 


2 


0 4 


2 


1 


0 4 1 


2 


0 4 


2 


1 


4 0 


1 


2 


4 1 


2 


0 


4 1 


2 


0 


0 1 


2 


3 


3 1 


0 


2 


3 0 2 


1 


3 0 


1 


2 


3 1 


2 


0 


3 0 


1 


2 


0 3 


1 


2 


(2) 






(1) 






CLSb 




CLSc 




(8) 






(7) 






(10) 






2 4 


0 


3 


2 4 


0 


3 


2 4 0 


3 


3 0 


4 


2 


3 0 


4 


2 


3 0 


4 


2 


3 0 


4 


2 


3 1 


4 


0 


3 1 


4 


0 


4 3 1 


0 


0 4 


1 


3 


4 1 


0 


3 


4 1 


0 


3 


4 3 


0 


1 


4 0 


1 


2 


4 0 


2 


1 


0 1 4 


2 


4 1 


2 


0 


0 4 


2 


1 


2 4 


1 


0 


2 4 


1 


0 


0 3 


2 


1 


0 3 


1 


2 


3 0 2 


1 


2 3 


0 


1 


2 3 


1 


0 


0 3 


2 


1 


0 1 


2 


3 


(6) 






(9) 






(4) 




(2) 






(6) 






(8) 






(3) 






3 0 


4 


2 


3 0 


4 


2 


3 4 0 


2 


3 4 


0 


2 


3 4 


0 


2 


3 4 


0 


2 


3 4 


2 


0 


4 3 


1 


0 


4 3 


1 


0 


0 1 4 


3 


0 3 


4 


1 


4 0 


1 


3 


4 3 


0 


1 


0 1 


4 


3 


0 4 


2 


1 


2 4 


0 


1 


4 0 2 


1 


4 1 


2 


0 


2 1 


4 


0 


2 0 


4 


1 


4 0 


1 


2 


2 1 


0 


3 


0 1 


2 


3 


2 3 1 


0 


2 0 


1 


3 


0 3 


2 


1 


0 1 


2 


3 


2 3 


0 


1 


(1) 






(5) 






(2) 




CLSd 




(4) 






CLSe 




(4) 






3 4 


2 


0 


3 4 


2 


0 


3 4 2 


0 


3 4 


2 


0 


3 4 


2 


0 


4 0 


2 


3 


4 0 


2 


3 


0 3 


4 


1 


0 3 


4 


1 


4 0 1 


3 


4 1 


0 


3 


4 3 


0 


1 


0 3 


4 


1 


3 1 


4 


0 


4 0 


1 


2 


4 1 


0 


2 


0 1 4 


2 


2 0 


4 


1 


0 1 


4 


2 


2 4 


1 


0 


0 4 


1 


2 


2 1 


0 


3 


2 0 


1 


3 


2 3 1 


0 


0 3 


1 


2 


2 0 


1 


3 


3 1 


0 


2 


2 3 


0 


1 


(3) 






(5) 






(2) 




(1) 






(8) 






(9) 






(10) 






4 0 


2 


3 


4 0 


2 


3 


4 0 2 


3 


4 0 


2 


3 


4 3 


0 


2 


4 3 


0 


2 


4 3 


0 


2 


3 1 


4 


0 


3 4 


0 


1 


3 4 0 


1 


3 4 


1 


0 


0 1 


4 


3 


0 4 


1 


3 


0 4 


1 


3 


2 4 


0 


1 


0 1 


4 


2 


2 1 4 


0 


0 1 


4 


2 


2 4 


1 


0 


2 0 


4 


1 


2 1 


4 


0 


0 3 


1 


2 


2 3 


1 


0 


0 3 1 


2 


2 3 


0 


1 


3 0 


2 


1 


3 1 


2 


0 


3 0 


2 


1 


(6) 






CLSf 




(7) 




(7) 






(6) 






(2) 






(3) 






4 3 


0 


2 


4 3 


0 


2 


4 3 0 


2 


4 3 


2 


0 


4 3 


2 


0 


4 3 


2 


0 


4 3 


2 


0 


3 0 


4 


1 


3 1 


4 


0 


3 4 1 


0 


0 1 


4 


3 


0 4 


1 


3 


3 U 


4 


1 


3 4 


0 


1 


2 4 


1 


0 


0 4 


2 


1 


2 0 4 


1 


2 4 


0 


1 


2 0 


4 


1 


0 4 


1 


2 


0 1 


4 


2 


0 1 


2 


3 


2 0 


1 


3 


0 1 2 


3 


3 0 


1 


2 


3 1 


0 


2 


2 1 


0 


3 


2 0 


1 


3 
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APPENDIX B 

NON-CLS PERMUTATIONS SATISFYING (13.2) 



For the permutations below, any integer y can be added to all terms, provided that reduction 
modulo k follows. Any term can be taken as the initial term, if followed by the k-1 terms to the 
right in order; thus each solution below implies k'^-l additional solutions. 



k = 7; 
k = 9; 



k = 11; 



k = 13: 



0 


2 


6 5 


3 


1 4 






















0 


2 


5 1 


6 


4 3 














0 


2 


5 7 


1 


3 8 


6 


4 


















0 


2 


6 1 


7 


4 8 


3 


5 










0 


3 


6 2 


1 


8 7 


5 


4 


















0 


4 


1 7 


6 


3 2 


8 


5 










0 


2 


7 5 


1 


8 4 


6 


3 


















0 


2 


4 8 


7 


3 1 


6 


5 










0 


3 


8 2 


7 


6 1 


5 


4 


















0 


2 


7 1 


6 


8 5 


4 


3 










0 


2 


5 7 


1 


4 8 


3 


6 


















0 


2 


6 1 


7 


2 5 


8 


4 










0 


6 


3 1 


8 


7 5 


4 


2 


















0 


4 


1 8 


6 


2 7 


5 


3 










0 


3 


6 1 


7 


4 2 


8 


5 


















0 


2 


5 8 


1 


7 4 


6 


3 










0 


3 


7 2 


8 


6 4 


1 


5 


















0 


3 


1 8 


7 


6 4 


2 


5 










0 


2 


7 


6 


1 


9 


5 


3 


10 


4 


8 










0 


2 


4 


9 


7 


3 


10 


1 


5 


8 


6 


0 


2 


9 


8 


6 


3 


10 


4 


7 


1 


5 










0 


3 


10 


9 


8 


4 


7 


5 


2 


1 


6 


0 


2 


10 


7 


6 


3 


9 


1 


4 


8 


5 










0 


2 


7 


1 


6 


9 


3 


10 


4 


8 


5 


0 


2 


6 


10 


9 


4 


8 


5 


3 


1 


7 










0 


4 


1 


8 


2 


9 


7 


3 


10 


6 


5 


0 


2 


5 


10 


3 


9 


8 


1 


6 


4 


7 










0 


2 


5 


1 


10 


9 


8 


6 


4 


3 


7 


0 


2 


8 


5 


1 


10 


9 


6 


4 


7 


3 










0 


2 


4 


7 


1 


10 


9 


6 


3 


5 


8 


0 


2 


10 


6 


9 


7 


4 


3 


1 


8 


5 










0 


2 


4 


9 


7 


10 


3 


5 


1 


8 


6 


0 


2 


6 


9 


7 


1 


5 


4 


10 


3 


8 










0 


3 


8 


2 


1 


9 


7 


10 


6 


5 


4 


0 


2 


8 


5 


9 


1 


4 


10 


7 


6 


3 










0 


2 


8 


5 


9 


3 


10 


4 


7 


1 


6 


0 


4 


9 


5 


3 


10 


7 


2 


1 


6 


8 










0 


4 


10 


9 


5 


1 


8 


6 


2 


7 


3 


0 


2 


8 


1 


6 


9 


3 


10 


7 


S 


4 










0 


2 


6 


10 


9 


7 


5 


4 


3 


1 


8 


0 


2 


6 


1 


8 


5 


3 


10 


9 


4 


7 










0 


2 


4 


7 


10 


1 


9 


6 


5 


3 


8 


0 


2 


9 


1 


12 


8 


10 


6 


4 


11 


7 


3 


5 


0 


2 


10 


12 


8 


4 


11 


9 


5 


7 


3 


1 


6 


0 


3 


8 


7 


12 


2 


11 


5 


4 


10 


9 


1 


6 


0 


3 


8 


7 


12 


2 


9 


1 


6 


5 


11 


)0 


4 


0 


3 


7 


4 


12 


2 


10 


1 


11 


8 


6 


9 


5 


0 


2 


6 


9 


1 


4 


11 


3 


10 


12 


8 


5 


7 


0 


3 


7 


11 


8 


1 


12 


5 


9 


6 


4 


10 


2 


0 


4 


9 


8 


2 


11 


5 


3 


12 


10 


7 


6 


1 


0 


4 


3 


11 


10 


7 


5 


1 


12 


6 


2 


9 


8 


0 
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APPENDIX C 

AN 18 X 18 GRECO-LATIN SQUARE 
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The above 18x18 square was constructed by the following strategy: 

(1) Obtain some odd number a, k/4 < a < k/3. 

(2) Partition k xk array S into rectangles Rj|^,R 2 f R^f R^: 

i) Rji^ = first k-a rows, last k-a columns of S. 

ii) R 2 = first k-a rows, first a columns of S. 

iii) R 3 = last a rows, last k-a columns of S. 

iv) R^ = last a rows, first a columns of S. 

(3) Let P(F) and P(f) be permutations F^, F 2 ,.. . fi,f 2 » . • . .^k-a 2(k-a) types alto- 

gether. 

(4) Let P(E) and P(e) be permutations Ej^ , E 2 , . . . , E^ and ej^, 02 , . . . , e^of 2 a types different from 
the types in (3). 

(5) Write F^ into A, where A is some randomly chosen set of k-a cells lying on k-a different 
rows, columns, and diagonals (as defined in (13. 1)) of R]^. 

( 6 ) Write P(F) into the diagonals of R]^ as in (13. 4). 

(7) Write P(E) into the unoccupied cells of the first row of R^. 

( 8 ) Write Ej^ into every cell of R]^ with Ej^ on the same diagonal by (7). 

(9) Write Fj^ into some cell of every row of R 2 not having an F^ in R^ by (5), so Fj lies in every 
column of R 2 . 

(10) Write Fj^ into some cell of every column of R^ not having an Fj^ in Rjj^ by (5), so F^ lies in 
every row of R^* 

(11) Write P(F) cyclically into the columns of R 2 and the rows of R 3 , from top to bottom, left to 
right. 

(12) Call the last k-2a rows, first k-2a columns of Rjj^ ”R*5'. " 

(13) Find z(c) setsC^ of k-2a cells in R’S* such that 

i) every cell in Cj^ lies on a different row, column, and diagonal of R]^; 

ii) every type in P(E) is contained once in 

iii) k-3a different types in P(F) occur in 

(14) Call the set of types Cj^DP(F) " 

(15) Let the first a rows and a columns of R 2 be called R>!‘=’5'. 

(16) Let the last a rows and columns of R 3 be called R=’!‘>!«*. 

(17) Find some set of 2a cells, a in R**!* and a in R**’!'; 

i) every cell in lies on a different row, column of S. 
ii) Bjj^ contains 2a different types from P(F). 
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(18) Call the types in P(F) not in 

(19) If for some Cj^, go on to (24). 

(20) If for n = 1,2, . . . ,z(c) go back to (17) and find B 2 # satisfying definition of 

If necessary, go on to B^, B^, . . . , B^. If = Y^, for some B^, C^.^, go on to (24). 

(21) If for m = 2,3,..., z(b) no = Y^^^, go back to step (5) and write F j into A 2 # A, where 

defined as A. If necessary, go on to » A^. If Xj^ = Y^^^ for some A^, B^^, Cj^, go 

on to (24). 

(22) If for u = 2,3,..., z(u) X^^ # Y^^, go back to (1) and obtain a 2 # a, a 2 defined as a. If neces- 
sary, go on to a 3 , a^, . . . , a^. If X^ = Y^ for some a^, A^, B^, C^, go on to (24). 

(23) If for V = 2, 3, ... , z(v) X^^ # Y^, give up. 

(24) If X = Y^, write f, into the cells of B^ and C_. 

' n ni I m n 

(25) Write P(f) cyclically into the diagonals of containing fp as in (6). 

( 26 ) Write P(e) into the cells of the first row of unoccupied by fj^, as in (7). 

(27) Write constantly into diagonals of Rj^ containing ej^, as in (8). 

(28) Write P(f) cyclically into the columns of R 2 and into the rows of R 3 as in (11). 

( 29 ) Construct by Algorithm I or II an a-order square having P(E) in the cells as types, ortho- 

gonal to a square having P(e) as types. Write both squares into the cells of R^. 

(30) Now S is a Greco-Latin square of order k". 

For the 18x18 square exhibited above, the solution was obtained when v = l,u = l,m = 4, z(c) 
= 12. For k" = 18, a is uniquely = 5. This solution took an hour by hand. Note that the values 
z(b), z(u), and z(v) are arbitrary limits, set as numbers of trials or as amounts of computer 
time. Of these four, z(v) will tend to have the smallest limiting values and z(u) the largest. In 
the opinion of this writer, it would be worthwhile to test z(v) and z(c) exhaustively, if possible, 
contenting oneself v'ith a fraction of the possible values for z(u) and z(b). 
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